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In the vicinity of the quark-hadron critical point, in the phase diagram of QCD, 
simple power-law relations constrain the mid-rapidity net-baryon density profile, 
for different heavy-ion processes, in a unifying scheme. The corresponding scaling 
variable specifies the proximity of a given experiment to the critical point. An 
important phcnomenological aspect of the critical properties in the baryonic sector 
is the possibility to observe intermittency effects in the phase space distribution of 
net baryons. 



1 Introduction 

Recent investigations concerning the QCD phase diagram (p — T) predict the 
existence of an endpoint along the critical line of the first order quark-hadron 
phase transition 111. This critical point, located on a line of nonzero baryonic 
density p — p c , is of second order and belongs to the universality class of a 3 — d 
Ising system. The QCD critical point communicates with the hadronic world 
through the fluctuations of a scalar field (cr-field) which carries the quantum 
numbers of an isoscalar (cr-meson) as the manifestation of a quark condensate, 
a ~ (qq), in thermal environment. In the effective theory of the QCD critical 
point, the classical cr-field is a natural order parameter, the fluctuations of 
which obey scaling laws dictated by the critical exponents of the 3 — d Ising 
system (r) rs 0, /3 R* 8 m 5, v rs |). In a baryonic environment, however, 
the chiral condensate is expected to have, at T = T c , a strong dependence on 
the net-baryon density, driving the cr-field close to zero for p m p c : (qq) p ~ 

^ ( P ~p Pa ) (Q9)o + 0[(p~ Pc) 2 ] where A is a dimensionless constant of the order of 
unity. This dependence suggests a new order parameter, m — p—p c , associated 
with the critical properties of the baryonic fluid created in a quark-hadron 
phase transition. In fact, approaching the critical point in the phase diagram, 
both the cr-field fluctuations and the fluctuations of the order parameter m(x) 
obey the same scaling laws ((qq) p ~ m(x)). 

In the next section we will exploit the scaling properties of the order param- 
eter m(x), properly adjusted to measurable quantities, in heavy ion collisions. 
For this purpose we consider net-baryon production in collisions of heavy nu- 
clei with total number of participants At and initial energy corresponding to 
a total size in rapidity Ay = L. 
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2 Scaling of the baryonic density 

The created baryons in the process of quark-hadron phase transition occupy a 

1/3 

cylindrical volume with transverse radius R± ~ and longitudinal size L 
(in rapidity). The parameter A±^ specifies the effective number of participants, 
contributing to the transverse geometry of the collision, and it is assumed 
Aj_ w 4f, valid both for central (A± = A rn i n ) and non-central collisions. 
Projecting out the net-baryon system onto the longitudinal direction we end 
up with a 1 — d liquid confined in a finite rapidity region of size L with local 
density p(y) — , directly related to the measurable net-baryon density 

in rapidity rife(y) = Putting — RqA 1 / 3 , we introduce a characteristic 
volume V Q — irR^Tf in terms of the freeze-out time scale Tt > 6 — 8 fm and the 
nuclear-size scale R Q which contains also any growth effects near the critical 
point (R > 1.2 fm). ). Using V~ x as a scale for baryonic, freeze-out densities, 
the order parameter m(y) of the 1 — d baryonic liquid is written: 

m(y) = Al 2/3 n b (y) - p c , 0<y<L (1) 

Near the critical point (T — > T c ) the order parameter m{y) obeys a scaling 
law: 

m(y) = iP [F (y/L) + tF 1 (y/L) + ...] (2) 

where t = T "~ T (T <T C ) and (3 is the appropriate critical exponent (/3 w 1/3). 
The leading term F D (y/L) in the series (||) is a universal scaling function 
whereas the nonleading terms Fi{y/L) [i = 1,2,..) are nonuniversal but A- 
independent quantities. Using eqs.(Q) and (|^) we may specify the freeze-out 
line in the phase diagram by considering the measurable, bulk density at midra- 
pidity rif, = n<b(L/2) as a function of the freeze-out temperature Tf. 

Al 2/3 n b = p c + t f j[F + t f F 1 + ...} (3) 

where tf = Tc ^ T/ and = Fj(l/2). Integrating now eq.(j|) in the interval 
< y < L we obtain at T = Tf. 

A^AL- 1 = Pc + t fi f {I +t f h + ...) (4) 

where Ii = Fi(£)d£. Introducing the variable z c = A ± 2 ^ 3 ' A t L~ Y we find 
from eqs.(Q) and (^) a scaling law for the net-baryon density at midrapidity: 

A 2l3 n h = ty(z c ,p c ) ; z c > p c (5) 
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where the scaling function ^>(z c ,p c ) has the property ^(z c — p Cl p c ) = p c . In 
the crossover regime z c < p c , where critical fluctuations disappear, the local 
density at midrapidity is, to a good approximation, rib ~ A t L~ x suggesting a 
continuous extension of the scaling law (||) in this region with ^(z c ,p c ) = z c 
(z c < p c ). It is of interest to note that although the scaling function ^>(z c , p c ) 
is continuous at the critical point (z c — p c ), the first derivative is expected 
to be discontinuous at this point in accordance with the nature of the phase 
transition (critical point of second order). Keeping only the next to leading 
term F\{y/L) in (j|) we finally obtain: 

A- A 2/3 n b = Pc + ^[f(z C7 p c )f + 

+C[f(z c ,p c )f +1 (6-a) 

A- ± 2/3 n b = p c + F o t f (l+t f ^—) (6-b) 



f(Zc,Pc) = + y/l + 2G(z c -p c y/P); 

z c > Pc (6-c) 



where C=^V,G=^. 

The component F a (y/L) which dominates the order parameter m(jj) in the 
limit T — > T c is approximately constant in the central region (y « , L ^> 1), 
far from the walls (at the points y = 0, L) due to the approximate translational 
invariance of the finite system in this region. On the other hand, approaching 
the walls, F a {y/L) describes the density correlation with the endpoints and 
obeys appropriate power laws. The solution which fulfils these requirements 
is: 

Fo{o= am -or 0/v (o<e<i) (?) 

and the final form of the scaling function ^(z c , p c ) becomes: 



40 /» 

*(z c , p c ) = p c + £ 1-2.) ^ <yZc ' Pc " 

\ V 1 V ' 



+C[f(z c ,Pc)f +1 ; z c >pc (8-a) 
^>(z c ,Pc) = z c ; z c < p c (8-b) 



where f(z c ,p c ) is given by eq.(6-c). It is straightforward to show that the 
discontinuity of the first derivative of \& (z c , p c ) at z c = p c is a nonzero universal 

constant: disc (gj-) ~ 1 — — (77 ~ §) as expected from the characteristic 

properties of a second-order phase transition. 
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3 Approaching the critical point in experiments with heavy ions 

Combining eqs.(^]-8) one may propose a framework for the treatment of cer- 
tain phenomenological aspects of the QCD critical point. The scaling law (|j) 
combined with eq.(8) involves three nonuniversal parameters: (a) the critical 
density p c and (b) the constants C, G which give a measure of the nonleading 
effects, allowing to accomodate in the scaling function (8) processes not very 
close to the critical point. We have used measurements at the SPS in order to 
fix these parameters on the basis of eqs.(||) and (8). More specifically, in a se- 
ries of experiments (Pb+Pb, S+Au, S+Ag, S+S) with central and noncentral 
(Pb+Pb) collisions at the SPS net baryons have been measured at midrapidity 
whereas the scaling variable z c = A ± 2 ^ 3 A t L~ x , associated with these experi- 
ments, covers a sufficiently wide range of values (1 < z c < 2) allowing for a 
best fit solution. The outcome of the fit is consistent with the choice G ~ 
and the equations (6) are simplified as follows: 

A~ 2/3 n b = p c + - (z c - p c ) + C(z c - p c ) 4 (9-a) 

4-2/3 . 2/ i/3 ttCI^ 3 
A x n b = p c -\ 1 (l + t f — - — ) (9-b) 



K\/y{L - y) 

+0[(z c - p c f] (9-c) 

In eqs.(9) we have used the approximate values of the critical exponents (3 « |, 
^ w i in the 3 — d Ising universality class. We have also added eq.(9-c) which 
gives the universal behaviour of the net-baryon density n b (y), in the vicinity of 
the critical point (z c — p c <C 1). The fitted values of the parameters in eq. ( |9-b[ ) 
are p c = 0.81, C — 0.68 and the overall behaviour of the solution is shown 
in Fig. 1. In turns out that the critical density is rather small, compared 
to the normal nuclear density p w 0.17 /m -3 : p c < ^? (R > 1.2 fm, 
T f > 6 fm), suggesting that the critical temperature remains close to the 
value T c a 140 MeV obtained in studies of QCD on the lattice at zero chemical 
potentially. 

The difference d c = \z c — p c \ in Fig. 1 is a measure of the proximity of 
a given experiment (A t ,L) to the critical point. We observe that the central 
Pb+Pb collisions at the SPS (d c f» 1.1) drive the system into the most distant 
freeze-out area, from the critical point, as compared to other processes at the 
same energy. In fact, the most suitable experiments to bring quark matter 
close to the critical point at the SPS are: S+S (d c rj 0.18), S;+Si (d c w 0.20) 
and C+C (d c as 0.06) ), central collisions. 
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4 Critical fluctuations in the baryonic sector 



Once the deconfined phase of quark matter has reached the critical point in 
a particular class of experiments, as discussed in the previous section, strong 
critical fluctuations are expected to form intermittency patterns both in the 
pion and net-baryon sector. As already mentioned in the introduction, the 
origin of these fluctuations can be traced in the presence, at T = T c , of a 
zero mass field with a classical profile (er-field) which, under the assumption 
of a phase transition in local thermal equilibrium, is described by an effective 
.action in 3 — d, the projection of which onto rapidity space is written as follows 



Ca 



dy 



da 
dy 



2c A py +i 



Ca = ^, Pc^T- 1 (10) 

Pc 



Equation (|10|) gives the free energy of the cr-field within a cluster of size 5y 
in rapidity and R± in transverse space. The critical fluctuations generated hy 
( |l2| ) in the pion sector have been studied extensively in our previous work El, 
therefore, in what follows, we are going to discuss the fluctuations induced by 
the cr-field in the net-baryon sector, noting that a direct measurement of these 
fluctuations becomes feasible in current and future heavy-ion experiments. For 
this purpose we introduce in eq.(|l^) the order parameter m(y) through the 
following equations: 



<r(y) « F(3 2 m(y) 



1 / dm 



r c ~ gi / dy 

ISy 



2 \8y 



g 2 \m 



s+i 



(qq)o ~ -3 fm 3 

m(y) = (3lm(y) (11) 



where: gi = F 2 g-i = 2C 2 A F 4 . The partition function Z = j V[m} e - r ^ 

for each cluster is saturated by instanton-like configurations! which for Sy < 5 C 
lead to self-similar structures, characterized by a pair-correlation function of 
the form: 



(rh(y)rh(O)} w — 



5T(l/3) fwR 2 ± C A 



6T(l/6) V P 2 



F- X y~ 



(12) 



2/3 



The size, in rapidity, of these fractal clusters is 5 C « ( 1 J^p2 J according 

to the geometrical description of the critical systems!. Integrating eq. (|l2|) we 
find the fluctuation (Snt,) of the net-baryon multiplicity with respect to the 
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critical occupation number within each cluster, as follows: 



is v ra F -i / ^jC^ 2 a / a r(l/3) 6/6 

<fab) ~ F V~W) r(i/6) V (13) 

The dimensionless parameter F 1 is of the order 10 2 and the size 5 C , on general 
grounds (R± ~ 2t c ) is of the order of one (<5 C ~ 1). As a result the global bary- 
onic system (in RHIC the size of the system is L ss 11) develops fluctuations 



at all scales in rapidity since the direct correlation (12) propagates along the 
entire system through the cooperation of many self-similar clusters of relatively 
small size (S c w 0.35 and (Sni,) « 140). We have quantified this mechanism 
in a Monte-Carlo simulation for the conditions of the experiments at RHIC 
(L ~ 11) in order to generate baryons with critical fluctuations. The distribu- 
tion of these "critical" baryons in the rapidity space for a typical event as well 
as the corresponding intermittency analysis in terms of factorial moments are 
presented in Fig. 2. 

The intermittency exponent of the second moment F2 in rapidity is found 
to be S2 ~ 0.18 which is very close to the theoretically expected value (g) of a 
monofractal 1 — d set with fractal dimension |. In conclusion, wc have shown 

6 " 

that measurements of net-baryon spectra in rapidity provide a valuable set of 
obsevables in heavy ion experiments, in connection with the phenomenology 
of the QCD critical point. 
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Figure 1: The scaling law (8) is illustrated together with measurements at the SPS. The 
critical point and the corresponding break in the slope of ty(z c ,p c ) are also shown. 
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(a) 




ln(M) 



Figure 2: (a) The fluctuations of critical baryons in rapidity for a MC-generated event, (b) 
The first three factorial moments for the event shown in (a) in a log-log plot. A linear fit 
determining the slope S2 (~ 0.18) of the second moment is also shown. 
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